Noise limits in matter-wave interferometry using degenerate quantum gases. 
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We analyze the phase resolution limit of a Mach-Zehnder atom interferometer whose input consists 
of degenerate quantum gases of either bosons or fermions. For degenerate gases, the number of atoms 
within one de Broglie wavelength is larger than unity, so that atom-atom interactions and quantum 
statistics are no longer negligible. We show that for equal atom numbers, the phase resolution 
achievable with fermions is noticeably better than for interacting bosons. 
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Atom interferometers jjj offer a potential sensitivity 
that exceeds that of their optical counterparts by as much 
as Mc 2 /Tiuj ~ 10 10 (|, p|. Besides being useful for funda- 
mental tests of quantum mechanics and precision mea- 
surements of fundamental constants, this property offers 
great promise for inertial navigation as well as for geo- 
physical applications, tests of general relativity, etc. j|, ||] 
fn particular, recent years have witnessed significant ex- 
perimental progress towards the development of atomic 
rotation sensors ^, Q with sensitivities now surpass- 
ing those of optical gyroscopes [0. So far, experiments 
with atom interferometers in the spatial domain have al- 
ways been carried out in the nondegenerate regime where 
quantum statistics and many-body interactions are neg- 
ligible and the interference pattern is the sum of single 
atom interference patterns. This is true even in those sit- 
uations where Bosc-Einstein condensates (BEC's) have 
been used as sources: in such cases, the interference 
experiments were performed after the atoms have been 
released from their trapping potentials and allowed to 
undergo a ballistic expansion in order to eliminate any 
mean-field effects, see for example jl0[ O. Likewise, 
the theoretical analysis of the sensitivity of atom inter- 
ferometers has also been confined to this regime. Refs. 
H ||, H analyzed the quantum noise limit of a matter- 
wave Mach-Zchndcr interferometer and showed that for 
uncorrelated inputs the phase resolution is 1/yN, for 
bosons as well as for fermions. Here, N is the total num- 
ber of particles entering the interferometer. The phase 
resolution limit results from the intrinsic quantum fluctu- 
ations associated with the sorting of the particles by the 
beam splitters in the interferometer. The reason that 
bosons and fermions yield the same sensitivity in these 
treatments is a consequence of the assumption that the 
fermions are quasi-monochromatic, i.e. the fermions all 
occupy the same state. This can only be justified when 
the atomic wave functions do not overlap [3|. 

Very much like the invention of the laser revolution- 
ized optics and interferometry, it would seem desirable 
to extend the operation of matter-wave interferometer 
to truly quantum-degenerate atom sources. It is there- 
fore important to extend these theoretical results to the 
case of highly degenerate beams of bosons or fermions 



for which p\d > I, where p is the atomic density and Xd 
is the de Broglie wavelength. In this limit, the effects of 
quantum statistics and of atom-atom interactions are no 
longer negligible, and it is expected that the sensitivity 
of the interferometer will decrease. For bosonic atoms, 
the interference pattern and phase sensitivity is degraded 
because of the atom-atom interactions, which give rise to 
a nonlinear phase shift as the atoms propagate along the 
arms of the interferometer. In contrast, for spin-polarized 
fermions — where s-wave scattering is absent for T ~ 
- the phase sensitivity is degraded because the Pauli 
exclusion principle forces each of the atoms in the inter- 
ferometer to occupy a different momentum state. Our 
central result is that contrary to what one might expect, 
the sensitivity of a Mach-Zehnder interferometer using 
a degenerate beam of fermions can be noticeably better 
than what can be achieved using a BEC. 

Figure 1 shows the atomic Mach-Zehnder interferom- 
eter that we consider, with input ports A and B. For 
convenience we assume that it can be treated as an effec- 
tive one-dimensional system in which the propagation of 
the atoms transverse to the interferometer arms can be 
ignored. That is, we assume that the atoms are tightly 
confined in the directions transverse to the interferom- 
eter paths, with the same wave function <f>o(r±). This 
holds when p < hu>x where p is the chemical potential 
for the bosons or fermions and htoj_ is the spacing be- 
tween the energy levels of the transverse wave functions, 
assumed to be harmonic. We introduce |he annihilation 
(creation) operators dk (dk ) and bk {bk ) for atoms en- 
tering the input ports A and B, respectively, with mo- 
mentum hk directed along either the upper or lower in- 
terferometer path. They obey the commutation relations 
a k a\, ± a\,a k = Sk,k> an d b k b\, ± b\,bk = 5k,k' where the 
minus sign is for bosons and the plus sign is for fermions. 

The beam splitters 1 and 2 are identical 50-50 beam 
splitters, which in the case of the first beam splitter has 
the action 
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where a Ui k(ai,k) are the operators for the upper (lower) 
paths of the interferometer. The mirrors Mi and Mi 
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FIG. 1: Mach-Zehnder interferometer. £ u and £i are the path 
lengths for the upper and lower arms, respectively. 

each result in a ir/2 phase shift. The Hamiltonian gov- 
erning the evolution of the Heisenberg operators in the 
two branches of the interferometer is 

s=l,u k ^ " ' 

where u> k = hk 2 /2M and M is the atomic mass. We re- 
call that x = for spin-polarized fcrmions at T ~ 0, while 
for bosons x = (Anha/Ml) J c£rj_|</>o(r_i_)| 4 , a being the 
s-wave scattering length and £ » ^; )U . The Hamiltonian 
ignores momentum changing collisions, which for an ini- 
tial state consisting of a weakly interacting BEC amounts 
to neglecting the quasiparticle fluctuations about the 
macroscopically occupied condensate mode. This is 
equivalent to a single-mode approximation for the con- 
densate. 

After the beams are recombined at the second beam 
splitter, the atoms are counted at detectors located at 
the output ports C and D. The relevant observable is 
the difference in the number of counts, 

M = N D -N C = Y^ (d\d k - c£c fc ) (f ) 

k 
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where d k and c k are the annihilation operators for the 
output ports D and C, respectively. 

In what follows, we consider only uncorrelated inputs 
where the initial quantum state can be written as the 
factorized product of Fock states \^a) \^b)- In this case, 
we recall that for nonintcracting atoms occupying the 
same state, the expectation value of N is (Af) = (Na — 
Nb) cos 9 where 9 is the phase difference between the two 
interferometer paths and Na and Nb are the number of 
atoms entering via the two input ports. Hence, N is a 
direct measure of 9, with fluctuations 

A9= AAA (2) 
\d(M)/d9\ 

and AA" = [(Af 2 ) - (Af) 2 } 1 / 2 §, |, g, We now com- 
pare A9 for the cases of a single mode BEC and of a 
spin-polarized Fermi gas. 

Bose-Einstein Condensate: Here we consider Na 
bosons entering port A and Nb entering port B in a time 
t, each atom having velocity vq = Hko/m. As such, (J\f) 
represents the number of atoms counted in the interval 
t. The quantum state describing this incident beam is 

l*(0)> = (4 ) NA (b{ ) NB \0)/VNaW.. 

We evaluate the moments of N by using an angu- 
lar momentum representation for the bosonic operators 
© @, 2J Z = a{ ( a ko - b{ o b ko , 2J X = a{ o b ko + ^a fco , 

2iJ y = d\b ko - b\ o a ko , and J 2 = N/2(N/2 + I) where 
N = a ko a ko + b\J) ka . In this representation, the incident 
beam state is |^(0)) = \j,m), where 

j = (N A + N B )/2, m=(N A - N B )/2. 

After some algebra one finds 



(AO = 2md^] m (P) cosfl + ( y/(j - m)(J + m + l)d^ +1 , m (/3) + Vti + ™)(j - m + l)$_ 1>m (/3)) sintf, (3) 

(A -2 ) - j(j + I ) + m 2 + f 3 . m (W) cos(20 + ft + g j<m (20) sin(20 + 0), (4) 

where 

f j>m (2(3) = (3m 2 - j(j + l))dg] m (20) - \^(j - m + 2)(j + m - I)(j - m + + ri)<$_ 2im {20) 

- \V(j +™ + 2)(j - m - l)(j + m + - m)d < £ +2 ^(2(3), (5) 

g j<m {20) = (2m - 1) V(i - m + l)(j + m)^ ) _ 1>m (2/3) + (2m + 1) V(j + m + l)(j - m)d% +hm (2(3), (6) 



and d^l m ((3) = (j,fn'\ exp(—if3j y )\j,m) are the matrix elements of a rotation about the j/-axis by an angle f3 = 
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x(Pl + *«) ■ I n these expression, is again the phase 
difference between the two interferometer paths, 

9 = {u ko + w_l + x(AT A + AT S - l)/2) {U-tl)+Pfa (7) 

where ij = 2£i/vq ( t M = 2£ u /vq) is the transit time 
for a point of constant atomic phase through the lower 
(upper) arm of the interferometer. Note that when \ = 

0, d~^, m (0) = 5 m >,m and our results agree with prior 
calculations p2| . 

The nonlinearity has a two-fold effect on the interfer- 
ence pattern. This is seen by noting that the phase ac- 
quired by an atom propagating along the upper or lower 
path is 4> s = (u ko +ui + Xa s>ko a St k )t s , and includes 
a self-phase modulation due to the other atoms travel- 
ing along the same path. However, the wave functions of 
the atoms entering the interferometer are coherently split 
between the two paths so that the effect of the nonlin- 
earity is no longer a simple self-phase modulation when 
expressed in terms of the incident quantum fields, 

<i>u~4>i = (w fco + oj± + xN/2)(t u - t t ) - x(*« + k)J y - 

The phase difference includes a self-phase modulation 
proportional to half the total number, a simple conse- 
quence of the fact that the atoms are equally split be- 
tween the two paths. However (f> a also includes a term 
proportional to J y . This results from the fact that the 
wave functions of the atoms entering the interferometer 
are coherently split between the two paths at the first 
beam splitter. Hence, the J y contribution is a result of 
the indistinguishability of the bosons that enter via the 
two input ports. 

We remark that 9 = ($ u - + (3/2 - x(t u - U)/2 
where the last two terms, (3/2 — x(tu ~ ti)/2, can be 
interpreted as a "vacuum effect", resulting from the non- 
commutativity of a s ,k and &1 kg . Its main impact on 9 is 
to replace N = Na + Nb with N — 1, showing explicitly 
that the nonlinear phase shift vanishes if only one atom 
is present. The (3/2 simply shifts the fringe pattern and 
is included in the definition of 9 for convenience. 

It is clear from Eqs. (|J) and ([I]) and the definition 
of d^, m {(3) that j3j is the critical parameter that deter- 
mines the effect of the atom-atom interactions on A9. 
Relevant values of (3 can be estimated by assuming a 
typical chemical potential of n/% ~ 10 3 — 10 4 s _1 and 
that 0o( r J_) is the harmonic oscillator ground state for 
w± ^ Taking the perimeter of the interferometer to 
be £i + £ u ~ 1 cm and typical recoil velocities, vo ~ 1-10 
cm/s, gives (3 > 10 -3 for 87 Rb or 23 Na, so that (3j > 1 
typically. 

The results of a numerical evaluation of Eqs. (|^) 
and (Q) are shown in Figs. 2 and 3. We find that 
A9 is minimum at 9 = (2ri + l)7r/2 and exhibits 
sharp resonances at 6 = nix where n is an integer. 
These resonances disappear in the limit f3 — > 0. 
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FIG. 2: AO at 6 = vr/2 as a function of (3 for N A = 1000 and 
iV fl = (solid line), N A = 1000 and N B = 50 (dashed line), 
N A = 1000 and N B = 100 (dash dot line), and N A = 950 and 
iVs = 50 (dotted line). 
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FIG. 3: AO at 9 = tt/2 as a function of 1/\/Wa for N B = 
and (5 = (solid line), (3 = 0.0005 (dashed line), j3 = 
0.001 (dashed dot line), j3 = 0.0025 (dotted line), (3 = 0.005 
(circles). 



These features agree qualitatively with a pertur- 
bation expansion for AO to lowest order in f3j, AO = 
VUU + 1) - ™ 2 ] / 2 ^ 2 I 1 + W) 2 (a(m, j) + b{m, j) cot 2 0)) 
where a(m,j) and 6(m, j) are positive definite and of 
order unity. Fig. 2 shows that for fixed A^^ the nor- 
malized phase fluctuations, A0((3)/ A9((3 — 0), increase 
with increasing Nb and for fixed N, the nonlinearity 
also results in an increase in the phase fluctuations 
with increasing Nb- Figure 3 shows that for (3 < 10 -3 , 
AO continues to scale as 1/\/Na when Nb = while 
for (3 > 10 -3 , A9 can actually show an increase with 
increasing Na- 

Fermi Gas: As with the bosons, we consider uniform 
beams of A^t and Nb fermions entering the interferom- 
eter via input ports A and B in a time r at average 
velocities vq. Their densities are ha = Na/vqt and 
Ub = Nb/vqt, and the quantum state for the incident 
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fermions is then given by 

i*(o)>= n 4 n 5 ii°> 

\k-k \<k FA \k-k \<k FB 



where kp a — ^n A and kp B — i^n B are the Fermi wave 
numbers for the two beams. We note from the definitions 
of kpA and kpB that there exists a new fundamental 
length scale, hp — vqt that is important for the descrip- 
tion of the fermions but was unimportant for the bosons. 
It represents the length of the two incident beams and 
can be interpreted as the quantization volume for the 
fermions. 

(TV) and AAf 2 are now found to be 

(AT) = cos9[n A S(n A 9/2) ~n B S(n B 9/2)]/a 
AN 2 = (n A - n B - cos26» [n A S{n A 8) - n B S(n B 8)]) /2a 

where S(x) = sin(x)/x, a = (2n/ko)/Lp, 9 — k (£ u — 1\) 
is the phase difference obtained for a single atom with 
momentum hkg, and we have taken N A > N B with- 
out loss of generality. We have also introduced the di- 
mensionless densities n A — 2 / nn A /ko = 2kp A /k^ and 
n B = 2nn B /ko = 2kp B /ko. They represent the number 
of atoms within the de Broglie wavelength Ao = (2w)/ko, 
so that quantum degeneracy corresponds to n > 1. Note 
that n is also a measure of the degree of monochromatic- 
ity of the incident Fermi beam. Earlier treatments of 
fcrmion interferometers dealt with quasi-monochromatic 
beams that correspond to n -C 1. 

We first restrict ourselves to the case n B = 0. The 
phase uncertainty is then 



V0 2 (l-cos20S(fi A fl)) 

X |cos 6 cos(6n A /2) - S(n A 6/2)[6 sin0 + cos 6]\ ,[ ' 

which shows that A9 ~ 1 / v 7 N A . However, Eq. (0) also 
indicates that A8 is an increasing function of 9 and hence 
of the path length difference. This is easily understood by 
noting that the individual phases contributing to the to- 
tal interference pattern are uniformly distributed over the 
interval [6-69,9+69], where 66 = k FA \£ u -£i\ = n A \9\/2. 
Consequently, for degenerate fermion beams with n A > 
1, the interference pattern is undetectable for large path 
differences, \9\ > 7r. This will continue to hold when 
n B ^Q. 

Figure 4 shows the minimum value of A9 for fixed n A 
and various n B plotted as a function of yfa ~ 1/y/N \ It 
is clear that the phase fluctuations continue to scale as 
1/vJV. Moreover, a comparison with Fig. 2 shows that 
increasing the number of atoms entering via input port B 
does not result in as dramatic an increase in A9 as with 
bosons. Whereas for bosons having only ~ 10% of the 
atoms enter via port B could result in an increase in Ad 
by a factor of 5 or more for moderate values of (3, having 
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FIG. 4: Minimum value of AO as a function of a 1 ^ 2 for ua = 
2.5 and ub = (solid line), ub ~ 1 (dashed line), ub = 1.25 
(dashed dot line), and fis = 1.5 (dotted line). 
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FIG. 5: Minimum value of A8 for Nb ~ and bosons with 
(3 = (solid line), bosons with (3 = 0.001 (dash dot line), 
fermions with ua = 1 (squares), fermions with ua = 2 (tri- 
angles), and fermions with tla = 3 (circles). 

more than a third of the atoms enter via B now results 
in an increase in A6* of less than a factor of 2. 

Fig. 5 provides a direct comparison of Ad for the BEC 
and Fermi beams when N B — 0. Remarkably, the phase 
resolution limit achievable with the Fermi gas is notice- 
ably better than what can be achieved with even a non- 
interacting BEC. Moreover, increasing the fermions den- 
sity enhances the phase resolution. Of course, this in- 
crease comes at the expense of a limited range of 9 over 
which the interference pattern is discernible. For BECs, 
by contrast, A9 is independent of 6 for (3 = and even 
for (3 =/= 0, A9 is nearly constant except in the vicinity of 
the resonances. 

In conclusion, we have examined the phase resolution 
limit for an atomic Mach-Zchndcr using degenerate Bose 
and Fermi gases. Our results indicate that one can attain 
better phase resolution using Fermi gases at the expense 
of a limited range of path differences over which there is 
a discernible interference pattern. 
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